Cosmic opacity: cosmological-model-independent tests and 
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Abstract 

With assumptions that the violation of the distance-duality relation entirely arises from non- 
conservation of the photon number and the absorption is frequency independent in the observed 
frequency range, we perform cosmological-model-independent tests for the cosmic opacity. The 
observational data include the largest Union2.1 type la supernova sample, which is taken for 
observed D^, and galaxy cluster samples compiled by De Filippis et al. and Bonamente et 
al., which are responsible for providing observed -Da- Two parameterizations, t(z) = 2ez and 
t(z) = (1 + z) 2e — 1 are adopted for the optical depth associated to the cosmic absorption. 
We find that, an almost transparent universe is favored by Filippis et al. sample but it is only 
marginally accommodated by Bonomente et al. samples at 95.4% confidence level (C. L.) (even 
at 99.7% C. L. when the r < 100 kpc-cut spherical /3 model is considered). Taking the possible 
cosmic absorption (in 68.3% C. L. range) constrained from the model-independent tests into 
consideration, we correct the distance modulus of SNe la and then use them to study their 
cosmological implications. The constraints on the ACDM show that a decelerating expanding 
universe with Q\ = is only allowed at 99.7% C. L. by observations when the Bonamente et al. 
sample is considered. Therefore, our analysis suggests that an accelerated cosmic expansion is 
still needed to account for the dimming of SNe and the standard cosmological scenario remains 
being supported by current observations. 

PACS numbers: 95.36. +x, 04.50.Kd, 98.80.-k 



I. INTRODUCTION 

The type la supernovae (SNe la) are observed to be fainter than expected from the 
luminosity-redshift relationship in a decelerating universe. This unanticipated dimming 
was first attributed to an accelerating expansion of the universe [l|, |2|. Although the 
existence of cosmic acceleration has been verified by several other observations, initially, 
there had been some debates on the interpretation of underlying physical mechanism 
for the observed SNe la dimming. For example, dust in the Milk Way and oscillation 
of photons propagating in extragalactic magnetic fields into very light axions had been 
proposed to account for the dimming [3, |j| . Any kind of photon number violation, such 
as absorption, scattering or axion-photon mixing, sensibly imprints its influence on the 
Tolman test J5|, which can be rewritten as a relationship among cosmological distance 
measurements known as the famous distance-duality (DD) relation p-l8j, 

^(l + z)- 2 = l, (1) 

where z is the redshift, D^ and D^ are the luminosity distance and the angular diameter 
distance (ADD) respectively. This reciprocity relation holds for general metric theories of 
gravity in any background, not just in that of the Friedmann-Lemaitre- Robertson- Walker 
background, and it is also valid for all cosmological models based on the Riemannian ge- 
ometry. That is, its validity depends neither on Einstein field equations nor on the nature 
of the matter-energy content. The DD relation plays an important role in modern cos- 
mology [9M12I]. and, in most cases, it has been applied, without any doubt, to analyze the 
cosmological observations. However, the reciprocity relation may be violated if photons 
do not travel on null geodesies or the universe is opaque. 

Fortunately, it is in principle possible to perform a valid check on the DD relation by 
means of astronomical observations. The basic idea is to search for observational candi- 
dates with known intrinsic luminosities as well as intrinsic sizes, and then determine their 
-Dl and -Da to test the Etherington relation directly. It is difficult for us to find objects of 
the same class with both known intrinsic luminosities and intrinsic sizes. Thus, a ACDM 
cosmological model is usually assumed when one performs tests by utilizing observed D^ 



or Da 13Ml6J] and the results show that there is no strong evidence for deviations from 



the standard DD relation. Recently, a model-independent method has been proposed to 



test DD relation by considering two different classes of objects, for examp 
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e, SNe la and 



201 ] . For a given 



galaxy clusters, from which D^ and Dp^ are determined separately 

ADD data, in order to obtain the corresponding D^ from SNe la, a selection criteria 

Az = |zciuster — ^SNe i a | < 0.005 is adopted. Using the phenomenological parameterized 

forms 

^(1 + z)- 2 = rj(z) , (2) 



with rj(z) = 1 + i] z and r](z) = 1 + yr-^cb and the data from Union2 SN la 21] and 
galaxy clusters, it was found that the DD relation can be accommodated at la C. L. for 
the elliptical (3 model [22J and at 3a C. L. for the spherical /3 model 23]. More recently, 
in order to avoid the bias from the redshift difference between D^ and ADD, methods, 
such as, the binning of SNe la 24], the interpolation 25[ and local regression [26[ from 
nearby SNe la points for a given galaxy cluster, were proposed and similar results were 
obtained. In addition, The DD relation tests by use of SNe la and gas mass function 
of galaxy clusters were carried out and similar results were also obtained 27H29J. So, 
overall, all the tests performed so far show that there is no strong indication of the DD 
relation violation. However, let us note that systematic uncertainties resulting from the 
morphological models of galaxy clusters and the redshift difference between D^ and ADD 
might exert influences on DD relation tests. 

If one considers that the photon traveling along null geodesic is unassailable, the DD 
relation violation most likely implies non-conservation of the photon number which has a 
mundane origin (scattering from dust or free electron) or an exotic origin (photon decay or 
photon mixing with other light states such as the dark energy, dilaton or axion J4J, ll4j). In 
this case, the flux received by the observer will be reduced and so the universe is opaque. 



If we assume that the flux from the source is decreased by a 
(observed) luminosity distance differs from the "true" one 



actor e t( - z \ then the inferred 
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where r is the opacity parameter which denotes the optical depth associated to the cosmic 
absorption. Initially, More et al. [32J studied the cosmic opacity by examining the differ- 



ence of the opacity parameter at redshifts z = 0.20 and z = 0.35, At = r(0.35) — t(0.2), 
where the difference of the observational luminosity distance (ADl i0 i, s ) at these two red- 
shifts was estimated from two subsamples of ESSENCE SN la 33j and the corresponding 
AZ^L.tmc was derived from the distance measurements of baryonic acoustic feature [34j in 
the context of ACDM. Assuming flat priors on Q\ and Qu in the range of < Q < 1, 
and uniformly spaced values of At G [0,0.5], they found that a transparent universe is 
favored (posterior probabilities of At peaked at 0) and At < 0.13 at 95% C L.. This 
method has been applied to investigate the homogeneity of the cosmic opacity in differ- 
ent redshift regions and the results suggest that the cosmic opacity oscillates between 
zero and non-zero values as redshift varies [35j, l36|. Later, Avgoustidis et al. [371, |38] 



carried out further studies by assuming an optical depth parameterization t(z) = 2ez or 
t(z) = (1 + z) 2e — 1 for small e and z < 1. There they took the standard luminosity 
distance in the spatially flat ACDM ((1 + z) 2 Da(z,Hm)) and the Union SN la [39J for 
-^L.true and -DL.obs respectively. In addition to the SNe la data, they also used the measure- 
ments of the cosmic expansion H(z) J40L|41[. By taking e G [—0.5,0.5], &m £ [0, 1] and 
Hq G [74.2 — 3 x 3.36, 74.2 + 3 x 3.36] |4l|] all uniformly spaced over the relevant intervals 
in a flat ACDM model and performing a full Bayesian likelihood analysis, they obtained 
a result e = -0.04+o;°7 (2a C. L.), which corresponds to an opacity At < 0.012 (95% 
C. L.) for the redshift range between 0.2 and 0.35, almost a factor of 11 stronger than 



the constraint obtained in Ref. 



32J. Recently, Lima et al. [42| reexamined this issue by 



confronting the luminosity distance which is dependent on two free parameters, i.e., the so 
called cosmic absorption parameter (a*) and the matter density (flu), with observations, 
using a subsample of Union2 SN la obtained by selecting SNe la with redshifts greater 
than cz = 7000km/s in order to avoid effects from Hubble bubble, and they found that 
the Einstein-de Sitter model (f2 M = 1) could be allowed at 68.3% (95.4%) C. L. in the 
case of a constant (epoch-dependent) absorption and concluded that a cosmic absorption 
may be responsible for the dimming of the distant SNe la without the need of an acceler- 
ated expansion of the universe. However, all these studies concerning the cosmic opacity 
assume a (flat) ACDM model, and are thus model-dependent. 

Here we propose another model-independent method to examine the cosmic opacity 



and investigate its possible implications for the cosmic evolution. If we assume that the 
violation of the DD relation is purely caused by the photon number non-conservation, 
then we can find out whether the universe is opaque by checking the possible violation 
of the DD relation. It should be emphasized that the cosmic absorption not only affects 
the luminosity distance measurements of SNe la observations as shown in Eq. ([3]), but 
also exerts influences on the angular diameter distance measurements determined from 
SZE+X-ray surface brightness observations 43l. 1441] . 



AT 2 



S- 
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(4) 



where ATcmb is the temperature change due to the Sunyaev-Zel'dovich effect (SZE) 
when the cosmic microwave background (CMB) photons pass through the hot intra- 
cluster medium and Sx is the X-ray surface brightness of galaxy clusters. The SZE 
spectra distortion of the CMB is determined by measuring the intensity decrements, A/, 
which is sensitive to the cosmic absorption. Additionally, the surveys of X-ray surface 
brightness are also sensitive to the opacity of the universe. Supposing the absorption is 
frequency independent in the observed frequency range (from microwave band to X-ray 
band), the "true" ADD connects the observed one measured in an opaque universe with 
-DA,truc — D A luster ■ e T . Thus, in actual calculations, the DD relation takes the following 



form: 

Now we will use the the 
galaxy cluster samples 
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argest Union2.1 SN la sample 1 |46| and the ADD data from 



23j to test, model- independently, the possible violation of the 



DD relation, which can be translated to a possible cosmic opacity. The observed D^ and 
-Da come from the latest Union2.1 SN la and galaxy clusters samples 22|, |23j, respectively. 
Actually, there are a number of inherent uncertainties in the selected astrophysical objects 
from which the observed Da are derived, e. g., the cluster asphericity [22] and the model 
for the cluster gas distribution [23]. In this paper, we consider the elliptical (3 model 
galaxy clusters sample [22J, spherical (3 model, r < 100 kpc-cut spherical (3 model and 



1 See Ref. U 



hydrostatic equilibrium model galaxy clusters samples |23| to investigate the impact of 
these inherent uncertainties on the cosmic opacity test. 

II. DATA AND CONSTRAINT RESULTS 

In order to place constraints on the cosmic opacity parameter r, we first parameterize 
it with two monotonically increasing functions of redshift, i.e., t(z) = 2ez and r(z) = 
(1 + z) 2e — 1 [37]] . These two parameterizations are basically similar for z <C 1 but they 
differ when z is not very small. As the data applied in our following analysis discretely 
distribute in the redshift range 0.023 < z < 0.890, our analysis we are going to perform 
may tell us something about the possible dependence of the test results on the parametric 
forms for r. To obtain t q i, s (z) determined by the following expression: 



Tobs(z) = - In 



Dl N 



(6) 



Df ustei (l + z) 2 

the data pairs of observed D^ and Da almost at the same redshift should be supplied. For 
the observed D^, the largest Union2.1 SN la is considered. Galaxy cluster samples, where 



the Da are obtained by combining the SZE+X-ray surface brightness measurements 43j, 



44j . are responsible for providing the observed Da- The first one, including a selection of 
7 clusters compiled by Mason et al. [47J] and a sample of 18 clusters collected by Reese 
et al. [48], was re-analyzed by Filippis et al. 22] by assuming an elliptical /3 model for 



the galaxy clusters. The second kind of samples are compiled by Bonamente et al. 23] 
with three different models for the cluster plasma and dark matter distribution, i.e., the 
spherical /3 model, spherical /3 model with r < 100 kpc cut, and hydrostatic equilibrium 
model. Therefore, the data derived from these three different models are adopted to 
check whether the cosmic opacity tests are sensitive to the model for the cluster gas 
distribution. The observed D^ are binned from the data points of Union2.1 SN la with 
their redshifts satisfying the certain criteria (Az max = kninster — ^sNe ia| max < 0.005) to 



match the observational data of the ADD samples 24j, |49_|. This binning method can 
minimize the statistical errors originating from the redshift difference between D^ and 
-Da- On the other hand, we alter Az max from 0.000 to 0.005 to ensure the number of the 
clusters that share the same SNe to be as few as possible, so as to reduce the dependence 



of opacity tests on the correlation of redshift-matched SNe. After obtaining t i, s (z) from 
these selected data pairs, we estimate the free parameters of a given parametric form by 
using the standard minimum \ 2 route: 

[r(z;p) -T ohs (z)] 2 



x {z;p) 



£ 



a 



(7) 



^obs 



where a Tobs is the error of r obs associated with the observed D L and D&, and p represents 
the free parameters to be constrained. The graphic representations and numerical results 
of the probability distribution of the opacity parameter e constrained from the model- 
independent tests are shown in Figures ([U [2]) and Table (JI|). These suggest that the 
dependence of test results on the above-chosen parameterizations for t(z) is relatively 
weak. Similar to the results obtained by examining the cosmic opacity in a particular 



redshift range (0.20-0.35) 32l. 1371 . |38[ and deforming the DD relation in terms of the 
cosmic absorption parameter 42J, we find, from Figure ([T]), that an almost transparent 
universe is also favored by the elliptical /3 model galaxy clusters sample 22J. For the 
ADD samples given by Bonamente et al. 23j, the results are shown in Figure [2j We 
find that the results are nearly not sensitive to the model of cluster gas distribution and 
a transparent universe can only be marginally accommodated at 95.4% C. L. (even at 
99.7% C. L. when the r < 100 kpc-cut spherical /3 model is considered). That is, all the 
constraints on the opacity parameter obtained from the Bonamente et al. samples prefer 
an opaque universe. These results are clearly different from what were obtained based on 



the ACDM in Refs. 



32 



42], where a transparent universe is obviously supported. In fact, 



these cosmic opacity test results are very similar to the previous model-independent tests 



for the DD relation 



IMS, 



24 



251 ] . However, our objective here is the cosmic opacity 



test with an assumption that the violation of the DD relation entirely originates from the 
non-conservation of photon number, rather than the DD relation test itself. 

In order to explore the implications of the cosmic opacity, let us transform the SNe la 
distance modulus in a transparent universe into that in an opaque one 



fht 



fx ohs (z) -2.5[loge]r{z) 



(8) 



and study the cosmological constraints resulting from this correction. Since the high 
redshift galaxy cluster data are absent in our discussion of the cosmic opacity, the distance- 
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modulus-modified SNe la used to investigate the cosmological constraints are cut down 
from Union2.1 with the criteria z < 0.784 and z < 0.890 when clusters in the Filippis et al. 
sample and the Bonamente et al. samples are applied, respectively. Since the dependence 
of cosmic opacity tests on the above-chosen parametric forms of t(z) is relatively weak, 
we only consider results obtained from the linear parametrization (r(z) = 2ez) in the 
following cosmological implication analysis. For the flat ACDM model, different from 
the methods used in Refs. 37|, |38|, |42|, we examine the probability distributions of Qm 



by considering the possible (68.3% C. L. range) opacity parameter e constrained from 
the previous model-independent tests. The results are shown in Figures ([31 HJ). We find 
that the opacity parameter e constrained from previous cosmological-model-independent 
tests impact slightly on the likelihood distributions of Qm and a universe with Q\ > 
is required to account for the dimming of the SNe la. This differs from the results in 
Ref. [42|, where the Einstein-de Sitter universe (Hm = 1) can be easily accommodated at 
68.3% and 95.4% C. L. for the constant and epoch-dependent absorptions, respectively. 

Without a spatially flat universe prior, we also investigate the ACDM with the cor- 
rected distance modulus of Union2.1 SN la by taking the opacity parameter e in 68.3% 
C. L. range constrained from the previous model-independent tests into consideration. 
The linear parametrization for the cosmic opacity is also considered. The results pro- 
jecting to the Qm — ^a plane are shown in Figure |5j Because of the "marginalization" 
of e, which somewhat weakens the constraints on parameters Qm and Q\, the statement 
that the expansion of universe is accelerating is less eloquent than the one in Ref. [I], 
which concludes that a currently accelerating universe is needed at 99.9% (3.9<r) C L. to 
agree with their distance measurements of SNe la. However, we find that a decelerating 
universe with Q\ = is only allowed at 99.7% C L. by observations when the spherical 
/3 model is taken into account. So, the standard cosmological scenario is still supported 
by observations, although current data may favor a universe with non-zero opacity, 



III. CONCLUSION AND DISCUSSION 

In this paper, by considering the luminosity distances provided by the largest Union2.1 
SN la sample together with the ADD given by galaxy clusters samples, we first exam- 
ine the possible cosmic opacity in a cosmological-model-independent way. Two redshift- 
dependent parametric expressions: t(z) = 2ez and t(z) = (1 + z) 2e — 1 are considered 
to describe the optical depth associated to the cosmic absorption. The results suggest 
that the tests of cosmic opacity are not significantly sensitive to the parametrization for 
t(z). For the ADD sample compiled by Filippis et al. 22] with an elliptical /3 model, 



we obtain that a universe with little opacity (almost transparent) is favored. For the 



ADD samples given by Bonamente et al. 23|, where three different cluster gas distri 



bution models are applied, the test results suggest that a transparent universe can only 
be marginally consistent with observations at 95.4% C. L. (even at 99.7% C. L. as the 
r < 100 kpc-cut spherical /3 model is applied). These are fairly different from the con- 



clusions in Refs. 



32 



42j . By considering the possible cosmic opacity (68.3% C. L. range) 
constrained from the previous model-independent tests, we obtain the corrected distance 
modulus of SNe la and then use them to investigate its cosmological implications. In the 
context of a flat ACDM, the likelihood functions of Qm are examined. The results are 
shown in Figures (0 H]) . We find that the opacity parameter e constrained from previous 
cosmological-model-independent tests have slight influences on the probability distribu- 
tions of Qm and a universe with Q\ > is required to account for the dimming of the SNe 
la. Discarding the condition of the spatial flatness, we display the corresponding plots 
in the Qm — ^a plane in Figure El We find that a decelerating expanding universe with 
Q\ = is only accommodated by observations at 99.7% C. L. when the spherical /3 model 
is considered. That is, a positive cosmological constant is still needed to account for the 
dimming of SNe la and the standard cosmological scenario remains being supported by 
current observations. 

Finally, it should be pointed out that, the presence of systematic uncertainties in obser- 
vations, especially ADD measurements using SZE+X-ray surface brightness observations, 
and the assumption of the frequency independency of absorption in the observed frequency 



10 



range in our analysis might result in some biases of our test results. Moreover, as for the 
DD relation test, the morphological models of galaxy cluster may also exert a remarkable 
influence on the tests for cosmic opacity. In fact, any conclusion that current data may 
favor a non-zero opacity should be backed up with a thorough analysis of these system- 
atics. Therefore, we may expect more vigorous and convincing constraints on the cosmic 
opacity within the coming years with more precise data, especially the ADD data, and 
a deeper understanding for the absorption in various wavelength bands and the intrinsic 
three dimensional shape of clusters of galaxies. 
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FIG. 1: Probability distribution functions of opacity parameter e obtained from the De Filippis 
et al. sample and Union2.1 SN la pairs for two parameterizations: t(z) = 2ez (blue solid) and 
t(z) = (1 + z) 2e — 1 (red dashing). 
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FIG. 2: Probability distribution functions of opacity parameter e obtained from the Bonamente 
et al. samples and Union2.1 SN la pairs for two parameterizations: t(z) = 2ez (blue solid) 
and t(z) = (1 + z) 2e — 1 (red dashing). The left, middle, and right panel represent results 
constrained from the isothermal f3, r < 100 kpc-cut isothermal (3 and hydrostatic equilibrium 
model respectively. 



gas distribution model 



t(z) = 2ez 



r(z) 



[\ + z 



,2e 



1 



Elliptical /3 model < - 0.000 

Spherical /3 model 
Spherical j3 model(r < 100 kpc-cut) 
hydrostatic equilibrium model e = 0.066 



+0.059+0.127+0.199 
-0.055-0.110-0.160 



n nsi+ - 046 +°- 100 + - 158 

u - uol -0.042-0.085-0.124 

n i 9n +0. 047+0.101+0.143 
u - lzu -0. 043-0.086-0. 114 



+0.037+0.079+0.123 
-0.035-0.070-0.102 



n ni/i+°- 071 +°- 145 + - 219 
u.uj.^_ ,069-0 .138-0 .203 



0.096 



0.140 

= 0.080; 



-0.058+0.114+0.169 
-0.056-0.107-0.154 

+0.054+0.110+0.148 



-0.052-0.102-0.142 

1-0.046+0.090+0.135 
-0.045-0.086-0.126 



TABLE I: Summary of the results for different optical depth parameterizations and cluster gas 
distribution models. 
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FIG. 3: Upper: The probability distributions of Qm in the context of flat ACDM when 
the absorptions model-independently constrained from the combination of De Filippis et al. 
sample and Union2.1 SN la are considered. The green, blue and red zonal regions represent the 
spans of S1m at 68.3%, 95.4% and 99.7% C. L. respectively. As the results are not sensitive to 
the parametric form of r(z), here the linear expression (t(z) = 2ez) is applied. Lower: The 
projections of the upper zonal regions in the Qm — e plane. The red dot (J7m=0.285, e=0.009) 
represents the best fit case. 
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FIG. 4: Upper: The probability distributions of Om in the context of flat ACDM when the 
absorptions model-independently constrained from the combination of Bonemente et al. sample 
and Union2.1 SN la are considered. The green, blue and red zonal regions represent the spans 
of VLyi at 68.3%, 95.4% and 99.7% C. L. respectively. As the results are not sensitive to the 
parametric form of t(z), here the linear expression (t(z) = 2ez) is applied. Lower: The 
projections of the upper zonal regions in the Qm — e plane. The red dot (Om=0.397, e=0.081) 
represents the best fit case. 
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FIG. 5: Marginalized regions at 68.3%, 95.4% and 99.7% C. L. in the Q M - ^A plane for the 
corrected data of subsamples of Union2.1 SN la by considering the observational constrained 
cosmic absorptions. As the results are not sensitive to the parametric form of t(z), here the 
linear expression (t(z) = 2ez) is applied. The red stars (Om = 1.0, Qa = 0.0) represent the 
Einstein-de Sitter universe. The left(right) panel corresponds to the results obtained from the 
De Filippis et al.( Bonamente et al.) sample. 
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